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Transport process in molecular chain in nonequilibrium stationary state is theoretically investi- 
gated. The molecule is interacting at its both ends with thermal baths which has different tem- 
peratures, while no dissipation mechanism is contained inside the molecular chain. We have first 
obtained the nonequilibrium stationary state outside the Hilbert space in terms of the complex spec- 
tral representation of Liouvillian. The nonequilibrium stationary state is obtained as an eigenstate 
of the Liouvillian which is constructed through the collision invariant of the kinetic equation. The 
eigenstate of the Liouvillian contains an information of spatial correlation between the molecular 
chain and the thermal baths. While energy flow in the nonequilibrium state which is due to the first 
order correlation can be described by Landauer formula, the particle current due to the second order 
correlation cannot be described by the Landauer formula. The present method provides a simple 
and perspective way to evaluate the energy transport of molecular chain under the nonequilibrium 
situation. 

PACS numbers: 05.60.Gg, 05.70.Ln, 44.10.+i 



I. INTRODUCTION 

Recently there have been paid much attention to 
nonequilibrium transport processes through molecular 
wire junctions. [1, 2 As the system is shrunk down to a 
smaller size than the mean free path of a carrier, such 
as an electron, a phonon, or an exciton, the ballistic 
quantum transport exhibits a characteristically different 
feature from the corresponding bulk properties. For ex- 
ample, it was discovered that the electronic or thermal 
conductance of sub-micron size one-dimensional chain is 
quantized at low temperature. [SHS When the system size 
is further shrunk down to the nanometer size, a dis- 
cretized energy level structure of a nanowire plays an 
important role in transport processes through the res- 
onance effect. Recent ultrafast nonlinear optical spec- 
troscopy has revealed the temporal behavior of relaxation 
process of the photo-excited molecular chain. [7',"F Also 
in biomolecules, such as alpha-helix protein and DNA 
strands, which consists of molecular chain structures, the 
transport process in nonequilibrium state have been ex- 
tensively studied in order to clarify the biofunctions of 
these molecules from a microscopic view point. j9HT3] 

In order to understand the transport processes in the 
nonequilibrium stationary state, we need to know first of 
all how to describe the nonequilibrium stationary state 
with which a physical quantity can be obtained as an 
expectation value of an observable. While it is obvious 
that a stationary state is described by a canonical dis- 
tribution in a thermal equilibrium, it is not so easy to 
find the explicit form of the stationary state for nonequi- 
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librium situation. Meanwhile a phenomenological Lan- 
dauer formula has been extensively used as a useful tool 
to evaluate a transport coefficient in a nonequilibrium 
stationary state without explicitly describing the station- 
ary state. There have been enormous works trying to 
make clear its microscopic foundation and its applicabil- 
ity in terms of nonequilibrium Green's function method, 
quantum Langevin method, and so on. |T| [T4142T] 

On the other hand, Prigogine and his coworkers have 
developed a theory of the nonequilibrium statistical me- 
chanics in terms of the complex spectral representation 
of the Liouville-von Neumann super-operator, or simply 
called LiouvilUan.\T2\ (23] They have clarified that the 
eigenvalue problem of the Liouvillian is classified into in- 
dependent subspaces according to the order of correla- 
tions based on the concept of the dynamics of the cor- 
relation where the fundamental object is a correlation. 
In a classical gas system, for example, each correlation is 
characterized by the dependence of the distribution func- 
tion on the wave number which characterizes the spatial 
correlation: The inhomogeneity is a correlation compo- 
nent which has a single non-vanishing wave number of a 
particle, while inter-particle correlation is a component 
which has several non- vanishing wave numbers of differ- 
ent particles. [24H26] 

This classification in terms of the correlation enables 
to transform the eigenvalue problem of the Liouvillian to 
the eigenvalue problem of the collision operator of the ki- 
netic equation in a correlation subspace. It is well known 
that the collision operator plays a central role in nonequi- 
librium statistical mechanics, as seen in Boltzmann equa- 
tion or Fokker-Planck equation. [22^, [27] A striking find- 
ing of the theory is that the spectrum of the collision 
operator is identical with that of the Liouvillian, which 
signifies the direct link between the microscopic dynam- 
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ics governed by the Liouvillian and the phenomenological 
kinetic theory. ^24l [26j |28] Since the cohision operator is 
non-Hermite operator, the eigenvalues can take complex 
values which reflects the dissipation of the system. The 
eigenstate of the Liouvillian is obtained from the eigen- 
state of the collision operator by operating a creation-of- 
correlation operator onto the eigenstate of the collision 
operator which represents the transition from a privi- 
leged correlation subspace to other non-privileged corre- 
lation subspaces.[24l[28] The eigenstate of the Liouvillian 
is then represented as a series of the correlation generated 
by the dynamics of correlation. 

There have been a few applications of the theory to real 
physical svstems. [23l [25l [29] but the transport processes 
in nonequilibrium situations have not been fully investi- 
gated yet in the context of dynamics of correlation. Our 
aim here is to apply the theory to the real mesoscopic 
system under nonequilibrium situation and to systemati- 
cally derive the transport quantities in a stationary state. 

In the present work, we consider a molecular chain cou- 
pled with different thermal baths at the both ends, where 
a quantum particle is conflned in the molecule and it 
ballistically transfers within the molecule. An exchange 
of energy occurs between the molecule and the thermal 
baths, while there is no particle exchange between them 
as shown in Fig{l] The nonequilibrium stationary state is 
obtained as the zero eigenstate of the Liouvillian which is 
represented by superposition of different order of correla- 
tions. A leading term is the non-correlation component, 
or called vacuum of correlation^ followed by the higher 
order correlation components. As will be shown, the in- 
homogeneity component of a single particle is created 
by the second order interaction from the vacuum of cor- 
relation. This is contrast to the case of a well-known 
gas system where the inhomogeneity component is dis- 
connected from the vacuum of correlation subspace by 
the interaction, thus is not created from the vacuum of 
correlation. [22] 

Once the nonequilibrium stationary state is obtained 
as a zero eigenstate of the Liouvillian, a physical quantity 
is evaluated as an expectation value of a corresponding 
observable in the stationary state. Since the nonequilib- 
rium stationary state is an eigenstate of the Liouvillian, 
the expression of a physical quantity is justifled from the 
microscopic dynamics without any phenomenological as- 
sumptions. 

We have derived an expression of the energy flow with 
use of the nonequilibrium stationary state, and found 
that the flrst order correlation is responsible for the en- 
ergy flow. We shall show that the expression of the 
energy flow is cast into the Landauer formula with a 
characteristic transmission function, because the flrst 
oder correlation is directly related to a collision operator 
which gives a transition probability between the molecu- 
lar states, represented by Fermi's golden rule. 

On the other hand, the physical quantity described 
by the higher order correlation cannot be reduced in the 
Landauer formula. As an example, we consider the parti- 



cle current which is induced by external thermal force, as 
in the mechanical force when a molecule is subjected to 
an external electric fleld. However, it is well-known that 
the treatment of the thermal force is in general much 
more complex than a mechanical force, because the ther- 
mal force stems from a many-body dissipative effect. [30] 
We shall reveal that the particle current is attributed to 
the second order correlation component. Since there is 
no direct connection between the second order correla- 
tion and the collision operator, we cannot cast it into the 
Landauer formula unlike the case of the energy flow in 
this problem. 

In Section |TT| we present a model Hamiltonian of the 
molecular chain coupled with different thermal baths at 
its both ends. The eigenvalue problem of the Liouvillian 



is presented in Section III to obtain the nonequilibrium 
stationary state as zero eigenstate of the Liouvillian. En- 
ergy flow is obtained with use of the nonequilibrium sta- 
tionary state and Landauer formula is derived in Section 
IV where we explain why the energy flow can be cast 
into the Landauer formula. In Section [V| the induced 
polarization, or its conjugate current, is evaluated as an 
example of the higher order correlation which cannot be 
reduced to Landauer formula. We illustrate in Section 
VT] a n application to DNA molecular chain. In Section 
VII| we give some concluding remarks. 

For readers who are not familiar with the complex 
spectral representation of the Liouvillian in Appendix [A] 
we introduce the Liouville space and show the explicit 
expression of the interaction in terms of the Liouville 
basis, respectively. The complex spectral representation 
of Liouvillian is summarized in Appendix [B] where the 
eigenstate of the Liouvillian is given as a functional of 
the eigenstate of the collision operator. With use of the 
explicit expression of the interaction, we derive the col- 
lision operator in Appendix [Cj In Appendix [P] we derive 
a formula which we use in Section IIIII 



II. MODEL 

We consider a one-dimensional molecular chain con- 
sisting in A/'-molecular units each of which contains a 
single bound state. This molecular chain is coupled with 
two thermal baths at the both ends. The total Hamilto- 
nian is written as 



H = Hm + Hb + qHmb 



(1) 



where Hm and Hb describe the molecular chain and the 
thermal baths, respectively, and the interaction is repre- 
sented by Hmb with a dimensionless coupling constant 
9- 

We consider that a quantum particle transfers in the 
molecular chain. The molecular Hamiltonian is then de- 
scribed by a one-dimensional tight binding Hamiltonian: 

N 

Hm = ^em\m){m\- ^ Jm,m'\rn){m'\, (2) 
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FIG. 1: One- dimensional molecular chain composed of A/" 
molecular units each of which possesses a bound state. The 
left- and right-end states have the energies of Si = £l and 
£^ — €r, respectively, while the other states have the same 
energy of £m = £o for m = 2, . . . , A/" — 1. A particle trans- 
fers between these states with the transfer integral of J. The 
molecular chain interacts with the two thermal baths with 
different temperatures of Tl and Tr at the both ends of the 
molecule. 



where £m is an energy of the bound state |m) at the m- 
th molecular unit as shown in Figjl] We denote the left 
end state as \L) = |1) and the right end state as \R) = 
\N). We assume that = eo in the middle of the chain 
for m = 2, . . . , A/" — 1, while the left and the right end 
states have different energies of sl and sr. The second 
term represents the particle transfer where we take into 
account only the nearest neighbor transfer; (m, m') in the 
second term denotes taking a sum of the nearest neighbor 
bound states. In the present work, we assume a constant 
value for the transfer integrals of Jm,m' = J foi" any m 
and m', and we take J = 1 as an energy unit which also 
becomes a unit of temperature. 

The eigenvalue problem of the molecular Hamiltonian 
Hm with the JV dimension is solved to obtain the JV 
eigenstates as 

Hm\Ej)=E-\Ej) (i = l,--- ,Ar), (3) 

where an eigenstate \Ej) is represented by 



AT 



1%) = Yl ^rnj\m) . 



(4) 



Hereafter in order to avoid a heavy notation, we simply 
describe \Ej) as 

Both of the left and right thermal baths are assumed 
to be a three-dimensional harmonic crystal described by 
a Debye model in a large box of volume . The Hamil- 
tonian of the thermal bath systems reads 



r=L,R q 



(5) 



where 6^^q (r = L^R) are the annihilation operators of 
the phonons of the thermal baths with the energy disper- 
sion given by cj^^q = c|q|. We take the box normaliza- 
tion with a periodic boundary condition for the thermal 
bath systems which gives the discrete wave vectors for 
the phonon mode as 



(6) 



where j is a three-dimensional integer vector and Aq = 
27t/Lb' In the large volume limit for the bath systems 
{Lb ^ oo), we have 



^Yl — ^ / * " ' ^^^'0 — ^ ^^^^ ' 



(7) 



where Q {l/Aqf. 

With use of the density of states of the thermal phonon 
system per volume Vph{uj)^ we can change the integral of 
q to the integral of uj as 



-y 



where Vph{uj) is given by 

( 47ru;^ 







for < < ujd 
otherwise . 



(8) 



In Eq.Q Debye frequency a;/) is given by 



(9) 



(10) 



where Nb is a number of the normal phonon modes. In 
the present work, we take to be large enough com- 
pared to the other parameters so that the value of ujjj 
will not affect the results. 

The unperturbed Hamiltonian Hq is then defined as a 
sum of Hm and Hb' 



Hn=H 



M ■ 



Hp 



(11) 



We consider the interaction between the molecule and 
the thermal baths represented by 



MB 



--^^«H.q|i?){i?|(&K.q + &kq) 



MB + ^Mb) 1 



(12) 



where Vr^^ = R) are the interaction potentials, and 
we assume that Vr,q is independent of Vt and q, i.e.. 



Vr- We consider that in Eq.(12) the molecular chain is 



interacting with a thermal bath only at the both ends of 
the chain so that the energies of the end states \L) and 
\R) are fluctuated by the interaction. With use of the 
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where \Y + rj/2) and {Y — rj/2\ are the eigenstates of 
Hm'- \j) = |y + vm and \f) = \Y - r]/2). Then \r], Y)) 
becomes an eigenstate of Cm' 



with the eigenvalue of 



(18) 



FIG. 2: Molecular level and the interaction scheme. The 
molecular eigenstates \j) are coupled with the thermal baths. 



E 



(19) 



Note the property 



eigenstates of Hm^ the interaction Hamiltonian Hmb is 
represented by 

r^L,R q 

- 

In FigM we draw the molecular eigenstates \j) of Hm 
with the energy Ej. These states are coupled with ther- 
mal baths through the coupling of the end state compo- 
nents. 



III. NONEQUILIBRIUM STATE AS AN 
EIGENSTATE OF LIOUVILLIAN 



The time evolution of the total system obeys the 
Liouville-von Neumann equation 



i^^pit) = Cp(t) , 



(14) 



where p{t) is a density matrix of the total system and C 
is the Liouvillian defined by Cp = [H^p]/h. We assume 
that in the initial state the left and right thermal baths 
are in thermal equilibrium with different temperatures 
Tl and Tr, respectively: 



exp [ - J2q fu^r,qbl^qbr,q/kBTr] 



{r = L,R) ,(15) 



where Zr = Ylq (l ~ ex.p[—huJr,q/kBTr]) ^ with ks the 
Boltzmann constant. Hereafter we take /c^ = 1. The 
thermal distribution p^^ gives the Planck's distribution 



nr{oj) 



1 



exp[-huj/Tr] - 1 



(16) 



In the present work, we shall use the Liouville space 
representation which is briefiy summarized in Appendix 
[A[ As a basis set for the particle system we shall introduce 
a Wigner representation defined by 



\V,Y))^\Y+^;Y 



2' 



){Y- 



2' 



(17) 



-r],Y 



^r],Y 



(20) 



By this definition, r] = and f] ^ represent the diagonal 
and the off-diagonal components of the density matrix of 
the particle, respectively. Note that the representation 
of 1 7^, Y)) corresponds to the Wigner basis representation 
of P)) of a Boltzmann gas system where k and P cor- 
respond to T] and y, respectively. 

Similarly the Wigner representation is defined for the 
phonon systems, we write the usual Wigner representa- 
tion for a (r, q) phonon mode as 

|z^r,q, A^r,q)) = |^r,q; <,q)) , (21) 

where 



_ _ / lyj _ ^r,q + ^r,q 

^r,q — ^r,q ^r,q ' -^^^,q — q 



The Wigner basis {N})) satisfies 



(22) 



Cb\M,{N})) = MM,{N})), (23) 



where {• • • } denotes a set of all the phonon normal modes 

and UUJ = Y.r=L,R Eq ^r,qC^r,q. 

The eigenstates of the unperturbed Liouvillian is 
then represented by a tensor products of \ri^Y)) and 
IM,™>as 



£o|r?,r))®|M,{iV})) = {^^,Y + vuj)\i^,YM{v},{N})) . 

(24) 

In terms of the Wigner basis of the eigenstates of Cq^ 
we can classify the Liouville space according to the or- 
der of correlations. For that purpose, we introduce the 
projection operators that specify the correlation compo- 
nents, such as a one-particle distribution of the particle 
P^^\ correlation between the particle and the phonon 
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systems and so on: 

^ ® |{0},{iv}))(({0},{iv}| , 

{N} Y 

(25a) 

{N} Y 

{0}^, iV,,q, {iV};))((l.,,q, {0}^, iV,,q, {iVj^l , 

(25b) 

{AT} y 

((^r,qi , ^r,q2 , {0}q, , ^'r,qi , ^'r,q2 , {^Iq^ ,q2 I . 

(25c) 

where {• • • }q^,q2,... means a set of all the normal modes 
other than qi , q2 , • • • . In order to consider the nonequi- 
librium stationary state which is obtained in the long 
time limit, it is appropriate to consider eigenvalue prob- 
lem of £, instead of solving Eq.(14) as an initial value 
problem. For the system coupled with the thermal baths 
with an infinite degrees of freedoms, C may have a com- 
plex spectrum. As shown in Appendix |Bj the complex 
eigenvalue problem of the Liouvillian is classified by the 
correlations: 

£|F(^))) = Zi'')|F(''))) , IF'^^^C = ((F(^)|4^) , (26) 



where (/i) is a combined index of (r^, {vY) in Eqs.(|25|) and 
a is an index of an eigenstate of a (/i)-subspace. The 

eigenstates of \F^^Jj and i!^F^^\ are the left- and right- 

eigenstates of C with the complex eigenvalues zi^^ for the 
total system composed of the particle and phonon sys- 
tems. In Eq.(26), we have used |-)) notation instead of | •)) 



for the total system consisting in the particle and phonon 
systems. We briefiy summarize the complex spectral rep- 
resentation of Liouvillian in Appendix [Bj and the reader 
could consult Refs. [23| l24| |26] for the details. Since the 
nonequilibrium stationary state is time independent and 
is achieved in the long time, we should seek for the zero 
eigenstate of C 

As have shown in AppendixjBj by acting the projection 
operators on Eq.(26), we have the eigenvalue problem of 



a collision operator given by 



where 



(27) 



(28) 



is a 'privileged component of and is a normal- 



r(M) 



ization constant which is given by Eq.(B17). 



Here, ^'^^^ is the collision operator familiar to nonequi- 
librium statistical mechanics. |22] [27l |30j [3j This oper- 
ator is associated to diagonal transitions between two 



states corresponding to the same projection operator 

(29) 



where Q^^) = 1 - V'^^\ In Eq.([29) 
1 



Q(/i)/:Q(M) 



(30) 



is called the creation-of-correlation operator, or simply 
the creation operator.\2A^ Eqs.(26) and (27) show that the 
spectra of C and ^ are identical. This directly links the 
microscopic dynamics to macroscopic kinetic processes. 

We can expand the creation operator in a series expan- 
sion of g as 



1 



OO / 

X5 [q^'^^CmbV^'^^ 



1 



Q(m)£oQ(p) 



5=0 



(31) 



The lower order expansion of the creation operator reads 



1 

1 



X - 



Q'-^'^CmbP'^^K (32a) 
Q^'''>CmbV^'''> . (32b) 



Substituting Eq.(31) into Eq.(29) we have the expan- 
sion of i!^^^\z) as 

*(A')(2) = p('')£op(A') + g^^/\z) + 0{g^) , (33) 

where 

^(^\z) = V^''^CmbQ^''^C[^\z) . (34) 

On the other hand, the non-privileged component 
(Q(^) -component) is obtained as a functional of the privi- 
leged component by using the creation operator as shown 



in Eq.(|B7|). The eigenstate IF^^^^J) is then written as 



|Ff ))) = ./Nf^(p('^^ + C^'^\Z^/^)) \uf)) . (35) 

Since the nonequilibrium stationary state |i^o^^)) ^ ^^^^ 
eigenstate of £, it is represented by the zero eigenstate 
of the collision operator \uq^-')) as 

\f!,'^)) = \/iVr + C('\+iO)) lui'-^)) , (36) 
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where the direction of the analytic continuation of C^^^ 
is indicated hy z = +iO, which is consistent with the 
fact that the approach to equihbrium is achieved in our 
future. In the weak couphng case considered here, the 
eigenstate |i^o^^)) represented by the expansion of the 
interaction up to the second order as 



(0)n 



where C^^^ and €2"^ are given in Eqs.(l32|. 



r(0) 



(i + ^cT^(+io) 



•^^Cf (+zO) 



(37) 



As shown in AppendixjBj the left eigenstate is similarly 
determined by the expansion of the interaction up to the 
second order as 



(0), 



N, 



(0) 



1 



(38) 

where f>^\z) and V^\z) are respectively the first and 
second order destruction-of-correlation operators which 
are determined by the series expansion of g from V^^^ [z) 



in Eq.(B15). Moreover, ((i)o^^| i s the left eigenstate of the 



collision operator shown in Eq.( B18|, an d the normaliza- 



tion constant N^^ is given by Eq.(|B17). 



In order to determine \u, 



we shall now solve the 



eigenvalue problem of the collision operator Eq.(27). 
Since the number of degrees of freedom of the thermal 
phonon system is infinitely larger than the number of 
molecular system A/", it can be shown that the V^^^- 
component of the phonon systems does not change in 
time and stay in their initial canonical distributions char- 
acterized by the initial temperatures Tl and Tr. There- 
fore we search for the eigenstate in the form of the tensor 
product of the density matrices of the molecular system 
and the phonon systems as 



\uf)) 



(39) 



where \p^l P^r)) ^ tensor product of the Liouville space 
vector of the left and right thermal phonon equilibrium 
distributions represented by Eq.JlsL In Eq.(39), 
is an eigenstate of the reduced coTTision operalor for the 
molecular system defined by 



(40) 



where Til^r stands for taking a partial trace of the 
thermal phonon systems, and therefore ^^^\z) is still 
a super-operator working on the molecular system. 

Up to the second order of the interaction, the reduced 
collision operator is represented by 



1 



iO+ - £0 



(41) 



There are eight diagrams for ^^2^^ in our model as shown 
in Appendix [C) 



In the present model, the matrix element of ^2^^ is 
explicitly represented in terms of the transition probabil- 
ities as 



3,r|tf |0,F)) 



r=L,Rr]>0 



{{0,Y\¥°^\0,Y + r^))=ikJ^^+\ 
((0,y|*f |0,r-r,))=^fc,^■^-^ 

where k^'^^^ is a transition probability from 



(t? > 0) to the |0,y)) given by Eqs.(C12) and (C14). 



(42a) 
(42b) 
(42c) 

0,Y±7j)) 



When we define the reduced density operator as 



m ^ tvl^r \p{t))) 



(43) 



it is found that for the weakly coupled system the f{t) 
obeys the following Pauli master equation of 



d 
dt 



E E_ {^"^""/^ 
(*)}, 



r]y^Or=L,R 
_l,Y,Y+r]AO) 

where we have defined as 

-^{{0,Y\f{t))). 



(44) 



(45) 



It has been known that this type of the master equation 
has a unique zero eigenstate and the other eigenstates 
have negative imaginary values of their eigenvalues. [32 - 

35 The zero value right eigenstate of ^2^^ is called colli- 
sion invariant which is represented by 



(46) 



Here we take the normalization condition for (py as 

^0y = l, (47) 



which sives from Eq.J39| that Tr = 1. Corre- 

spondingly the collision mvariant for the left eigenstate 
is given by 



EE 

Y {N} 



which satisfies {{vj^^ k^'' ^ 



3,r|(({0},{7V}|, 



1. 



(48) 



Substituting Eq.(46) into Eqs.(39) and (37), the 



nonequilibrium stationary state is obtained up to the sec- 
ond order of the interaction as 



= v^(i + 5cr(+^o) + /cf(+^o)) 



Y r=L,R 



(49) 
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Similarly the left eigenstate for the stationary state can 
be obtained as 



an 



(0), 



Y {N} 



X (l + gV['^ (+iO) + g'vi'^ (+iO)) , (50) 
which satisfies {{F^^^\F^^^)) = 1 and 



Tr 



,(0) 



(0) 



(51) 



First we shall consider the case for A/" = 2 where the 
nonequilibrium stationary state is analytically obtained. 
The nonequilibrium stationary population (F = 1, 2) 
are obtained by 



n = 



^i(n^(A) + l)+4(n^(A) + l) 



^i(2n^(A) 



vlriLiA) 



■vli2nR{A) 



vl{2nL{A) + 1) + vl{2nR{A) + 1) 



(52a) 



(52b) 



where A is given as an energy differnce between the two 
molecular states of |£^2) and \Ei) (See Eq.(Il9|): 



E-2-Ei 1 



h 



h 



\/(eL - srY + 4J2 . (53) 



For a longer molecular chain, it is difficult to analyti- 
cally obtain the nonequilibrium stationary state. Instead 
we have numerically solved the eigenvalue problem of the 
collision operator and obtained the collision invariant. As 
an example, we show in FigjS] the results for a molecule 
with A/" = 10: the molecular level structure, the eigen- 
states of the molecular Hamiltonian (a), the spectrum of 
the collision operator (b), and the population of the col- 
lision invariant (c). In Figjsj^b) and (c) we have fixed at 
Tl = 2.0 and changed Tr and have taken vl = vn = 1 
and g = 0.1. Since there are ten basis states belonging 
to P(o) subspace ( |0,r)) with F = 1 to 10 ), we have 
ten eigenstates of the collision operator. For any Tr, 
there always exists an eigenstate with Zq^^ = 0, i.e. a 
collision invariant, while other eigenstates are decaying 
states with ImZ*^^^^^ < 0, as mentioned above. 

c), we have also shown by the dotted lines 
the canonical distribution for the average temperature 
Tm = {Tl-\-Tr)/2. At the low temperature, the nonequi- 
librium state population is different from the canoni- 
cal distribution for the average temperature, while both 
agree well at the higher temperature. Furthermore, at 
the low temperature case where the discrete molecular 
level structure plays a key role, the nonequilibrium state 
population cannot be represented by the canonical dis- 
tribution for any unique temperature. This suggests that 
due to the quantum effect the local temperature of the 
molecule cannot be identified at the low temperature 
case, in contrast to the classical system. [6J 
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FIG. 3: (a) The molecular level structure for the chain with 
J\f — 10, where the parameters are el = = 0, J = 1.0, vr = 
vl — V — 1.0, and g — 0.1. (b) The spectrum of the collision 
operator z'^^ as a function of Tr changing from to 4.0 with 
Tl — 2.0 fixed, where the vertical axis denotes — ImZi°\ (c) 
The state population of the nonequilibrium stationary state 
(j)Y as a function of Tr. The dotted lines are the population 
for the average temperature of Tm — {Tl + Tr)/2. 



IV. ENERGY FLOW 

In this section, we consider the energy fiow in the 
nonequilibrium stationary state in terms of the zero 
eigenstate of jC obtained in Section |III[ The energy fiow 
through the molecule is evaluated by considering the en- 
ergy change of the molecule. The energy change of the 
molecule in the nonequilibrium stationary state is ob- 
tained by 



w{t) 



d 
dt 



M)t 



dt 



Tr [HMp{t)] 



Using Eq.(|l4[), w{t) is represented by 

W{t) = -i{{HM\TTL<^R[Cp{t)])) 

= -ig{{HM\TTL^R[CMBp{t)])) 



(54) 



(55a) 
(55b) 



where we have used the fact that {CmHm)^ — and 
T^^L^nl^Bpit)] = in the second equality. Taking into 
account Cmb = ^mb + ^mb corresponding to Hmb 



TJfL 
^MB 



^MB Eq.(12), we can divide w{t) into two 
contributions due to the interactions with the left and 
the right thermal baths: 



wit) = w^{t)^w^{t) , 



(56) 
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where 

w'-it) = -ig{{HM\TrL^R [^mbPW] )) (r = L,R). 

(57) 

The sign of w^{t) is positive when energy flow comes into 
the molecule from the thermal bath r = L or R. 

Since >CFq^^ = in the nonequilibrium stationary 

state, we find that w{oo) = —i{{HM\T^^L^R [jCFq^^'\ )) = 0, 
resulting in w^{oo) = —w^{oo) which guarantees that in- 
fiux and outfiow of molecular energy are balanced in the 
stationary state. We then define the energy fiow from 
the left to the right thermal bathes going through the 
molecule as 



Since in the long time Eq.(51) indicates that 



|p(t)))^^^|Fo) 



N, 



(0) 



(58) 



(59) 



^y,y+77^^^^^ which is the transition rate per time be- 
tween the states of |0, Y)) and |0, F + r])) due to the cou- 
pling with the left thermal bath. The energy fiow coming 
in from the left thermal bath is obtained by multiplying 
with it the energy difference between the states of |0, Y)) 
and 10,1^ + 7^)) states, This correspondence is 

originated in the fact that the collision operator i^p^\z) 
can be represented by the first order creation operator 
C[^\z) as shown in Eq.(34). Therefore in the weak cou- 
pling case the energy fiow can be derived by using the 
solution of the eigenvalue problem of the collision op- 
erator: Simply multiplying the energy difference of the 
states and the transition rate between them. 

This correspondence naturally leads it to a Landauer 
formula which has been widely used to interpret the 
carrier fiow of the mesoscopic system in nonequilibrium 
situation. ^ Landauer formula reads 



hujT{uj){nL{uj) — nR{uj))duj , (62) 



it is found from Eq.(49) that U{oo) can be represented 
up to the second order of Cmb as 



U{oo) = -ig^{{HM\TTL^R 



(60) 

It should be noted that the first order correlation created 
from the collision invariant contributes to the energy fiow 
in the nonequilibrium stationary state. [22l |27l |3TJ Fur- 
thermore, by comparison with Eq.(34), it is found that 



the energy fiow is related to the collision operator which is 
expressed by the transition probability given in Eq.(42). 



Inserting Eqs.(46) and (32i) into (60) and using the 
matrix representation of given in Eqs.(A12), we 



can explicitly represent U^^ =U{oo) as a sum of any pair 
of the states of |0, Y)) and |0, F + //)): 

^st y^y^^UY^r],Y 

r]>0 Y 



r]>0 Y 



:|(/)ynL(A^^y+^) - 0y+^(nL(A^^y+^) + l)| 



(61a) 



2 5:^nA,,y+.{^W^y - fcr+>y+,} 



r]>0 Y 



(61b) 



where 2 ^^^q in Eq.(6l3) can be replaced by "^^^q- 
Here we would like to make a comment that the same 
formula can be derived by evaluating energy change of 
the thermal bath as shown in Appendix [P] 

We note that there is a striking correspondence be- 
tween Eq.([6l]3) and Eq.(44): The right hand side of the 
two equations involve a common factor of (/c^^^'^(/)y — 



where T{oj) is a transmission function determined by us- 
ing Fermi's golden rule. In the simplest approximation, 
T{cjo) = 1 is assumed, which results in a quantization of 
thermal conductance. [6 In more elaborate works, they 
have estimated the uj dependence of T{oj) refiecting a 
resonance effect due to the discretized molecular level 
structure. 17, 18 



In the present model, the equation (61 ) can be cast into 



the form of Eq.(62) in terms of the transmission function 
defined by 

r(c.)=2^^5(c.-A,,v+n) 

rj>0 Y 



nL{Ar,,Y+r,) " nii(A^,y+^) 



(63) 



The transmission function has a strong resonance char- 
acteristic for the molecular level structure, which is re- 
fiected in the delta function of Eq.(63): uj = {Ey+r] — 

EY)/h. 

While the energy fiow due to the first order correlation 
can be represented by the Landauer formula as shown 
above, we shall show in the next section a physical quan- 
tity due to the higher order correlation, such as an in- 
duced polarization, cannot be represented by the Lan- 
dauer formula. It is worthwhile to note that the Lan- 
dauer formula is derived in our approach through the 
resonance effects between the molecular chain and the 
thermal baths, i.e., the dissipation occurs at the edges of 
the system contacting to the baths. 

Before going to the next section, we shall show some 
examples of the energy fiow in the nonequilibrium sta- 
tionary state of the molecule. The energy fiow is analyt- 
ically obtained for a molecule with the length oi M = 2 
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FIG. 4: The energy flow U^^ of a molecule with the length of 
A/" = 10 in the Tl-Tr plane. The parameters for the calcu- 
lation are the same as used in Fig|3] The temperature range 
is taken from to 10 in (a), and from to 1 in (b). The 
energy flow as a function of Tr for several fixed values of Tl 
are shown in (c) and (d). 




" -^-^ " -3 -2 -1 



FIG. 5: The dependence of the energy flow U^* on sl and sr 
with Tl = 2.0 and Tr = 0.01 flxed, where the other param- 
eters are the same as used in Fig|4] The bird's eye view and 
the contour plot are shown in (a) and (b), respectively. 



U^^ is maximized at £l — ^i? = as a function of sl—Sr, 
because a transfer between the end states \L) and \R) 
through the molecule occurs most effectively. For differ- 
ent values of El and Sr^ the energy transfer is allowed 
to occur due to the energetically spread molecular states 
as shown in Fig|2j As temperature increases, the energy 
flow increases. 



by using Eq.(52). In this case, we have 



2 .,2 



Vr V 



(2nL(A) + l)vl + (2n^(A) + l)vl 



2 ' 



(64) 



yielding 



(2nL(A) + l)vl + (2nfl(A) + l)v% 



(65) 



where A is given by Eq.(53). 



For a longer molecule we have numerically calculated 
the energy flow. As an example, we have shown in Figj4] 
the energy flow through the same molecule as studied in 
Section III We show the energy flow U^^ in the bird's-eye 
view in the Tl-Tr plane in FigQa) and (b). The wide 
temperature region of U^^ is shown in (a), and the low 
temperature region in an expanded scale in (b). We also 
show the energy flow U^* as a function of Tr for various 
values of Tl in (c) and (d). It should be noted that while 
the energy flow is linearly proportional to the tempera- 
ture difference near thermal equilibrium Tl — Tr^ this 
linear relation breaks down in far-from equilibrium. 

The energy flow also depends on £l and £r. We 
show in Fig js] the energy flow as a function oi Sl and 
£r for = 2.0 and Tr = 0.01, where the other parame- 
ters are the same as in Figj4j The bird's eye view and the 
contour plot are shown in (a) and (b), respectively. The 



V. INDUCED POLARIZATION, PARTICLE 
CURRENT 

In spite of the fact that the thermal force stems from 
complicated many-body dissipative effect, when the po- 
larizable molecule is subject to an external thermal force 
under nonequilibrium condition, a polarization is in- 
duced, just like a simple mechanical force with an ex- 
ternal electric field. We shall show that an induced po- 
larization within the molecular states is attributed to the 
second order correlation so that the quantity cannot be 
reduced in a Landauer formula unlike the energy flow. 

Polarization operator is represented by 



ex 



\m){m\ 



(66) 



when Xm = md with a lattice constant d and the electric 
charge e. In the present work, we take e = 1. The 
conjugate current to the polarization is defined as a time 
derivative of the polarization: 



d 
dt 



Tr[xp{t)] . 



(67) 



Using Eq.([T4]), X{t) is then represented as 



-i{{x\TrL<^R[Cp{t)])) 
-i{{x\CM\m)) , (68) 



where we have used the fact that {{x\Ttl(^r[Cb + 
g^Ms] = 0, and f{t) is a reduced density operator of 
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the particle given by Eq.(43). When we define the cur- 
rent operator as 



X = i[HM,x] 



(69) 



then the current X{t) is given by X{t) = Tr[X/(t)]. Note 
that the current operator X may be represented by using 
the site basis {\m)} as 



N 



X = - ^ [i^M, \m){m\] 



id 



m=l 
N-1 



^ ^m,m+i (1^ + - \m){m + 1|) , 



(70) 

which agrees with the ordinary definition of the particle 
current in a one-dimensional discrete lattice. [36] 

Considering that pit) in the long time limit is given by 



F^^^ Eq.([59|, we then have 

X^' ^ -i{{x\TrL^R[CF^'\)) ={) , 



(71) 



i.e., the total particle current in the nonequilibrium sta- 
tionary state vanishes because CF^^^ = 0. 

Substituting the spectral representation of Cm given 
by Eq.(IT8l) into Eq.(|68l), we have 



At) = E PI^,nA,,y((r?,r|/(t))) 

(72a) 

Y r]>0 

x{{{n,Y\-{{-v,Y\)\f{t))) (72b) 
^EE^'^.yW. (72c) 

Y r]>0 



where we have defined a particle current component 
X^,y(t). The current component Xj^^yit) is a product of 
the three factors: i) Transition dipole moment between 
\Y + 77) and \Y — rf) molecular states, (lx\ri^Y)) = {Y — 
iri/2\x\Y ^ iri/2)^ ii) Transition frequency between these 
two states, A^^y, and iii) Off diagonal matrix element of 
the reduced density matrix, (((77,Y'| — (( — ri^Y\)\f(t))) . 
The product of the factors i) and ii) represents that the 
particle current is attributed to the transition of the par- 
ticle between \Y — iri/2) and |F + 7^/2) molecular states. 
Furthermore, since the factor iii) is written as 



(((r,,y|-((-r,,y| 
= {Y+l\f{t)\Y- 



\m)) 

l>-(^-ll/wll'+?>, 



(73) 



the reduced density matrix f{t) should be non-Hermitian 
in order to obtain a non- vanishing X^^y(t). This is 
quite contrast to a thermal equilibrium at /3 = 1/T 
where the density matrix is given by a Hermitian ma- 
trix Peg = ex.p[—/3H] which leads to a vanishing particle 
current component, while Z^V ^^^^ finite values in the 
nonequilibrium stationary state. 

Substituting Fq^^ into Eqs.(72 ), the current component 
is expressed by 



q-st 

V^Y 



= -i({i|r?,F))A,,y({{r?,F|-((-r?,F|) 

X 52|TVi«fi[cf (+iO)^(°V;i]», for > (74) 



It should be noted that the contribution to the current 
component X^^y is attributed to the second order creation 
operator. 



Inserting Eqs.(46) and ([32| into Eq.(74), the explicit 
expression of Z^V obtained as 



q-st 
r],Y 



or by using 



X E j ^Ph{^^,Y-r^-dW,Y-^-i?{^Y-^ (^r ( A^,y _^_^) + l) - 0y _ ^ _^n^ ( A^,y _^_^) } 

+^p/i(A^,Y-r7+^)|Cr,y-^+^P{^y-^^r(A^,y_^+^) - _ 5 +^ (n^ ( A^,y _^+^) + l) } 
+^p/i(A^,Y+r7-^)|Cr,y+^-^P{</>y+^ K ( A^,y +^_^) + l) - (/)y + ^ _^n^ ( A^,y+^_^) } 
+^p/i(A^,Y+r7+<e)kr,y+5+<eP{^Y+J^r(A^,y+^+^) - ^y+^+^ (n^(A^,y+^+^) + l) } 

the transition probabilities of given by Eqs.(C12), X^^y can be written as 



_R-kr (PY+^+0 



-kr ^ ^y+5+d 



(75) 



(76) 
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(a) y-r}l2 



Tl = 2.0 , Tr = 0.01 




FIG. 6: The particle current components X^^^q^ ^ for J\f ■■ 



10, 



VL = VR ^ 1.0, g 0.1, Tl 2.0, and Tr = 0.01. In (a), 
all the components of T^^^y ^tre drawn, and in (b) and (c) the 
components of 77 = 1 and 77 = 3 are shown, respectively. 
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FIG. 7: The temperature dependence of some current com- 
ponents shown in Fig|6] The current components as a func- 
tion of Tr for a fixed value of Tl = 2.0 are shown. The 
1^,Y for (7y,y) = (1,5/2), (1,9/2), (1,13/2), (1,17/2) and 
(77, Y) = (3, 7/2), (3, 9/2), (3, 13/2), (3, 17/2) are shown in (a) 
and (b), respectively. 



As seen in Eq.(|76|), in the calculation of the current com- 
ponent P^^^ molecular states be- 
tween \Y -\- r]/2) and \Y — r]/2) we have to consider the 
transitions between the other molecular state \Y + r]') 
and \Y) expressed by or , and between 

and \Y^r]) expressed by or 
As a result, the particle current cannot be cast into the 
Landauer formula contrary to the energy flow, where 
there is no other molecular states involved in the cal- 
culation of the energy flow l^Y-\-r],Y- 

As an example of the particle current we show in 
Figjo] the particle current components X^^y same 




FIG. 8: Space correlation in the molecule. All the parameters 
are the same as in Fig|6] (a) All the components of /m,M, and 
(b) the dependence on the relative distance m for M = 5 and 
M = 5.5. 



molecule which has been investigated in Section [rV] 
where we have taken JV = 10, vl = vr = 1.0, ^ = 0.1, 
Tl = 2.0, and Tr = 0.01. In (a), all the components of 
X^^Y are drawn, and in (b) and (c) the components for 
fixed values of 77 = 1 and 77 = 3 are shown, respectively, 
where the horizontal axis is Y. When £l = ~ ^ 

for an even integer of 77, because ((77, Y\x)) =0 due to the 
symmetry. It is found that Z^V takes a large value when 
T] = 1 which suggests that the particle current is large for 
a pair of adjacent molecular states in energy, and the X^^y 
becomes small as r] increases. It is also found that X^^y 
is positive for a large Y while it is negative for a small 
Y: The induced polarization due to the molecular states 
with a higher energy is directed from high temperature 
side to the low temperature side. 

We also show the temperature dependence of Z^V 
of FigjT] as a function of Tr with a fixed value 
of Tl = 2.0. In (a) are shown X^\ for 
(77, r) = (1,5/2), (1,9/2), (1,13/2), (1,17/2), and in (b) 
for (77, F) = (3, 7/2), (3, 9/2), (3, 13/2), (3, 17/2). The 
particle current vanishes at Tr = Tl, and they linearly 
depends on the temperature difference around Tr = Tl, 
while they change nonlinearly in the low temperature 
region reflecting the discrete molecular level structures. 
This behavior corresponds to that of the energy flow 
shown in Figj4j because both the energy flow and the 
particle current are born out from the unique collision 
invariant by the actions of C^^^ and C2^^ on |iio^^)). 

So far we have investigated the particle current com- 
ponents attributed to a quantum correlation between a 
pair of molecular states. With use of the representa- 
tion in terms of the site basis given in Eq.Q, we may 
reveal the quantum correlation in space which is gener- 
ated in the nonequilibrium stationary state. Similarly to 
Eq.([73|, we shall investigate the non-symmetricity of the 
off-diagonal elements of f{t) in terms of the site basis as 

/™,M = {{{m, M\-{{- m, M\) \f{t))) , (77) 

where \m^M)) is the Wigner basis in terms of the site 
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basis defined similarly to Eq.(17) by 

I ±m,M)) = \M± 



-;MTy) = |M±-)(MT 



(double sign in same order) 



2 ' ' 
(78) 



with use of site basis \M ± m/2). We show in Fig Js] the 
imaginary part of fm,M for the same case of Figj6j In (a) 
we show the all the components of fm,M and in (b) the 
dependence on the relative distance m for M = 5 and 
M = 5.5. As seen from the figures, the values of fm,M 
does not depend on the central position of M but on the 
relative distance of m, suggesting that the quantum cor- 
relation in space decreases, as m increases. We also found 
that /m,M = dit Tl = Tr: The quantum correlation in 
space comes to appear only at nonequilibrium situation. 



(a) 




(b) 




FIG. 9: The energy flow (a) and the particle current (b) of 
ome-dimensional base pairs of DNA molecule with 10 bases. 
In (b) Tl = 0.07 (350K) and Tr = 0.06 (300K) are used. 



VI. APPLICATION TO ONE-DIMENSIONAL 
DNA CHAIN 

As an application of our result, let us give an exam- 
ple of a real physical system, a hole transfer in one- 
dimensional chain of DNA base pairs. [9] In this case, a 
hole (or an electron) is an energy carrier and the molecu- 
lar system has been described by the tight binding Hamil- 
tonian given by Eq.([2|. It has been known that in a chain 
of DNA bases, the site energy of a guanine (G) and cy- 
tosine (G) base pair is higher than that of a thymine 
(T) and adenine (A) pair: sgc — ^ta — 0.5eV, and the 
transfer between the base pairs have been estimated to 
he J ^ 0.4eV. We show in Figj9]^a) a typical result of 
the energy fiow when we put two GG pairs at the both 
ends of a molecular chain with its length of A/" = 10, 
where we have taken ^ = 0.5. It is found that the energy 
fiow increases nonlinear ly with temperature; Note that 
Tl,r = 300K corresponds to T^ ji = 0.06 in our unit. 
The hole current in the molecular states is also shown 
in Figj9jb) and the schematic picture of the energy fiow 
and the particle current is shown in Figj9jc), where we 
fix Tl = 0.07 and Tr = 0.06 corresponding to 350K and 
300K, respectively. It is found that the energy is trans- 
ported by a particle current in the lower molecular states. 

VII. CONCLUDING REMARKS 

We have obtained the nonequilibrium stationary state 
under thermal force as the zero eigenstate of the Liouvil- 
lian of a molecular chain which is weakly coupled with 
different thermal baths at the both ends. The zero eigen- 
state is represented in the expansion of the order of cor- 
relation following the principle of the dynamics of corre- 
lation. The physical quantity in the nonequilibrium sta- 
tionary state is derived by taking an expectation value of 
an observable with respect to the stationary state. The 
energy fiow and the particle current are attributed to the 
first order and second order correlations created from the 
vacuum of correlation, respectively. Since the first order 



correlation is directly related to the collision operator in 
kinetic theory, the energy fiow can be cast into the form 
of Landauer equation with a transmission function with a 
strong resonance which refiects a discrete level structure 
of the molecular states. It is also found that the thermal 
force induces an polarization, or its conjugate particle 
current, which cannot be cast into the Landauer type 
formula, because the particle current is attributed to the 
second order correlation. Even so, there is a correspon- 
dence between them in their temperature dependences. 

Here we have dealt with the nonequilibrium transport 
process of a molecule coupling with a thermal phonon 
bath with very broad spectral width. It is interesting 
to investigate how the energy fiow will be changed if we 
modify the phonon density of states in such a way that a 
particular molecular states be resonantly excited. By re- 
placing a phonon field with a radiation field in the present 
work, we can investigate how the energy fiow behaves un- 
der a monochromatic light excitation in a nonequilibrium 
stationary state. These further extensions of the present 
model will be shown elsewhere. 
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Appendix A: Representation in Liouville space and 
the expression of Cmb 

In this section, we shall briefiy review the Liouville 
space representation of a Hilbert space operator. The 
Liouville space is spanned by linear operators in yl, 5, • • • 
in the ordinary wave function space. |24j As usual, the 
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inner product of the Liouville space is defined by 

{{A\B)) = TriA^B) , (Al) 

where A and B are linear operators acting on wave func- 
tions, and is a Hermite conjugate of A. For the case 
where the wave function space is spanned by a complete 
orthonormal basis, 



Y^\a){a\ = l,{a\(3) = 5, 



(A2) 



where |n^^q) is the number state for a thermal bath nor- 
mal mode (r, q) of Hb- In order to clarify the order of 
correlation, we prefer to use the Wigner basis defined by 
Eqs.(IT7l) and (pll): 



l^,n = IJi;j2)) 



where 



(A7) 



(A8) 



the Liouville space is spanned by a complete orthonormal 
basis of the dyads |a; P)) = |q^)(/3| , i.e., 

^|a;/3))((a;/3| = l, ((a; /^O) = (5,,a'(^/3,/3- (A3) 

The matrix element of the usual operator A in the wave 
function space is given by 



{{a;M = {a\A\f}) . 



(A4) 



The Liouville basis is constructed of the tensor product 
of the eigenstates of the unperturbed Hamiltonian. For 
the molecular system and the thermal bath systems, the 
Liouville basis are written by 



b'l; j2)) = \ji){j2\ , 



(A5) 



where \j) denotes the eigenstate of Hm given by Eq.Q, 
and 



l^r,q;<q)) = |^r,q)«q| 



(A6) 



and 



|^^r,q,A^r,q)) = |^r,q; <,q)) , 



(A9) 



where z/r,q and A^r,q are defined in Eq.(21). The Wigner 
basis is the eigenstate of Cq as shown in Eq.(24). The 



Wigner basis then form the complete orthonormal basis 
satisfying 

{(r?,F|r?',F')>® m.{m\{v'}.{N'])) = S,,,,6y,y, 



{u},{,.'}S{N},{N'} , 

E 



r],Y 



\V,Y))UY 



(AlO) 

m,{N}))m;{N}\ = i, 



(All) 



where IM, {AT})) ^n.,ql'^r,q,iV.,q)). 

Now we consider the matrix element of the interaction 
Liouvillian Cm b = Hm b x 1 — 1 x Hm b in terms of these 
Wigner basis. The calculation of the matrix elements of 
Cmb can be done in a straightforward manner, yielding 



ii^,Y\i{v}-{N}\gCMB\n',Y'))\W},{N'})) = j^ E E^'- 



{y + l\r){r\Y' + \)5y,_,^^y_jN,,^ + ^ + 1 --^•-^ ^K,^,^r,.+Au},W}hNUN'} 



r^L,R q 

1 d_ 



-i^' - j\r){r\Y - ^)V+v,y+zzAK,q - ^ -t,. ^^^^^^^^ j^^^^^^,^^^^^^^^,^ 



+ (r + ||r)(r|r + |)^^,_^^^_ W7V.,q + ^ 



V' 



+ 1 e^'^-^^<^,..,^-i^{,},{,/}^{iv},{iV'} 



(A12a) 
(A12b) 
(A12c) 
(A12d) 



where S' stands for the product of the Kronecker delta ex- 
cept for the interaction normal mode of (r, q). We draw 
the diagram of the correlation corresponding to these four 
terms in FigjTol In the figures, the dotted line denotes 
the correlation, i^^^q, of the phonon mode involved in the 
interaction, and the solid or double solid lines denote the 
correlation, 77 or 77', of the molecular state. The filled cir- 



cle stands for the vertex of the interaction whose matrix 
element is written at the vertex. The diagrams of (a) 
and (c) correspond to Eqs.(A12i) and (AI23), each of 
which is attributed to the 6^^q (phonon absorption) and 
6j! q (phonon emission) terms in Hmb x 1 term of Cmb^ 
where the phonon line appearing on the left side of the 
vertex implies that the transition of the particle ket state 
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^^,(y + ||r)(r|y + 7^'-|) 



(b) 1^'^)) l^r-,q,Ar,,q)) 



-9Vr{Y-v'^j\r){r\Y-j) 



(c) 



FIG. 10: The diagram of the interaction LiouviUian Lmb, where the Unes represent the correlation. [221 [24] The dotted hne 
denotes the correlation, z/r,q, of the phonon mode involved in the interaction, and the solid or double solid lines denote the 
correlation, 77 or 77^, of the molecular state. The filled circle stands for the vertex of the interaction whose matrix element is 
written at the vertex. 



of \Y' + 77Y2) to \Y -\-r]/2) occurs by the interaction. On 
the oth er ha nd, the diagr ams of (b) and (d) correspond 
to Eqs.( |A12| 3) and ( |A12| i), each of which is attributed 
to the q (phonon emission) and 6r,q (phonon absorp- 
tion) terms in 1 x Hmb term of jCmb^ where the phonon 



line appearing on the right side of the vertex implies that 
the transition of the particle bra state of {Y' — r]' /2\ to 
(Y — 77/2 1 occurs by the interaction. The reader may be 
referred to the textbook about the way of writing the 
correlation lines. ^22^ 



Appendix B: The Complex Spectral Representation 
of Liouvillean and Subdynamics 

In this section, we shall summarize the complex spec- 
tral representation of LiouviUian. [24 Useful formula for 
this paper are listed without proof. The reader may refer 
to some references for detail. j24H26] 

In the complex spectral representation of LiouviUian, 
we consider the eigenvalue problem for each correlation 
subspace (/i) = (77, z^), where 77 and v represent the order 
of the correlation of the particle and phonon, respectively. 

The complex eigenvalue problem is written as 



(Bl) 



where the LiouviUian can have complex eigenvalues 
ImZj^^ 7^ 0- It has been shown that the time evolu- 
tion splits into two semigroups; one is oriented toward 
our future t > with ImZ-^^ < (equilibrium is ap- 



proached for t 00), while the other is oriented toward 
our past t < with ImZj^^ > 0. All irreversible processes 
have the same time orientation. To be self-consistent we 
choose the semigroup oriented toward our future, which 
determines the direction of the analytic continuation of 
the eigenfunction of £.[ 24] [28] 

Now we introduce the projection operators defined in 



( 25 ) which satisfy 



We also introduce the projection operators 
which are orthogonal to 



(B2a) 
(B2b) 
(B2c) 



(B3) 
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We solve the eigenvalue problem (Bl ) for the perturbed Assuming completeness in the space V^^\ we may al- 
ways construct a set of states {(('^^•^^|} biorthogonal to 



system with ^ 7^ under the boundary conditions for the 
unperturbed case: 



{14"^))}, i-e., 



= hr^,,,- ^Y.\^fmf\ = 1 . (B12) 



(B4) 



Hence, Q^^)|Fj^^)) = for = 0. The Q^^) components 
are created through the interaction for ^ 0. The right 
and left eigenstates, and are biorthonor- 

mal sets satisfying 

^ ,^ = 1 . (B5) 



Applying the projection operators V^^^ and QS^^ in 



The equation (B9) combined with ( |B1Q ) shows that 



(25) and (B3) to the (Bl), we derive the set of equations: 



the V^^^ component of li^^ )) (which is called ''''privileged 

componenf of is an eigenstates of the collision 

operator, which has the same eigenvalue zj^^ as the Li- 
ouvillean. The solution of the eigenvalue problem of the 
Liouvillean for our class of singular functions has unique 
features. The privileged components satisfy closed equa- 
tions and the Q^^^ components are "driven" by the priv- 
ileged components [See (B7)]. 

Combining (B7) with (28), we obtain the right eigen- 



states of the Liouvillean given by 



Equation ( B6 3) leads to 



(B6a) 
(B6b) 



(B7) 



1^^ 



where 



))) ^ +C(^)(zj^)))|^x5^))) . (B13) 

Similarly, we obtain for the left eigenstates given by 
((Fj^^l = {{v^f\(^V^^^ ^V^^\Z^/^^ (B14) 

where the operator V^^\Z^^^) is called the destruction- 
of-correlation operator, or the destruction operator for 
short, and is defined by [cf. ([3Q|)] 



1 



1 



Q(/i)/:Q(M) 



Q^^^ . (B15) 



is called the creation-of-correlation operator, or simply 
the creation operator\2A^ Substituting (B7) into (B6i), 
we obtain 



where 



¥>'\z\^^)\uf)) = zf\uf) 



(B9) 



(BIO) 



and N^^^ is a normalization constant which will be deter- Eq.(B5)with Eqs.(B13) and (B14) as 
mined later. Here, ^^^^ is the collision operator familiar 
to nonequilibrium statistical mechanics. [22j [271 1301 El 
This operator is associated to diagonal transitions be- 
tween two states corresponding to the same projection 
operator P^^): 



Again V^^^ {z) corresponds to the off-diagonal transitions; 
f,{t^)(^z) = p(/^)p(M)(^)g(M). Using X>(^)(z), the collision 
operator ^^^^(z) is also written as 

+ V^^^V^^\z)Q^^^gCMBV^^^ . (B16) 
The normalization constant N^^^ is determined by 



+ V^^^gCMBQ^^^C^''\z)V^^^ . (Bll) 



Note that ( |B9| ) is a nonlinear equation in the same sense 
of the Brillouin-Wigner perturbation method, i.e., the 
eigenvalue Z-^^ appears in the collision operator. 



(B17) 

The states ((i)]^^ | are the left eigenstates of the collision 
operator 

= {{v\^\zf . (BIS) 

We have revealed the correspondence between the eigen- 
value problems of the Liouvillean C and the collision op- 
erator ^f^^\ 
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Appendix C: Collision operator in the vacuum of correlation subspace 



In this section, we shah give the exphcit expression of the cohision operator of ^2^^ by applying Eqs.(|A12|) into 



Eq.(41). We show the diagrams of ^2^\+i0) in Figjll] The diagrams (a), (b), (c), and (d) represent the loss of the 
stateJO, Y)) which are written by 

(a): ((0,y|*f |0,r))(,) = fji^ Yl E\''r\'KyCr,Y-,\ %.^ ^ — — K^.,) + 1) , (CI) 

(a') : ((0,y|*f |0,r))(,,) = fj^E E T.\''r\'KYCr,Y+,\ \.^ ^ —^UrM , (C2) 

(6): ((0,y|ff |0,F))(,) = p^E E El"-l'l<i-W.P ,n+ ^ -— n.KJ, (C3) 

(6'): ((0,y|*f |0,F))(,,) = fJi^ ^ Y.\^r\'KYCr,Y-,\ %.^ ^ -^KK.,) + 1). (C4) 



where nr(w) is Planck's distribution function given by Eq.([l6|. 

On the other hand, the diagrams (c), (c'), (d), and (d') represent the gain of the \0,Y)) state from the |0,y ± 77)) 



io,n 



io,n 









/ \-v,Y-vm 


A 

1 
t 

+ lr,q \ 


|-r/,r + r//2)) 


r ^ 


r 



I -77,y + r,/2)>| ;. 



|-77,r-,7/2» 



|o,y» 




io,n 




|o,r + r?» 




\o,Y-7^)) 


(a) 


(b) 


(c) 


id) 



|o,y» 



\v,Y + v/2)) 



|o,y» 



t 
t 

t 

\n,Y-T,/2))\/ 



\0,Y)) 
(a') 



\0,Y)) 



(b') 



io,n 







i 




\ri,Y-7j/2)) 





io,n 















\0,Y-v)) 



\0,Y + T,)) 



(C) 



(d') 



FIG. 11: Diagram of the collision operator ^2 (+^0) 
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states which are written as 



(c): 

(c') 

id): 



r|*(°)|o,y + 7?)) 



^ ^ |Ur|2|c*yCr,y+^| 



g r=L,R 



1 



1 



q r=L,R 



q r=L,R 



Z0+ - (A_^,y_^/2 +CJ^,g) 
1 

Z0+ - - UJr^q) 



{nr{0Jr,q) ^ I) ,(C5) 
^r(^r,g) , (C6) 
nr{uJr,q) , (C7) 
+ ,(C8) 



where n^(co') is the Planck's distribution given by Eq.( 16). The summation for the thermal phonon modes q is replaced 
by the integral with use of the density of states given by Eq.([9|. 

When we sum up the contributions of these diagrams, we have for > 



{{0,r|tf |0,F)) = -2Tri^ Y.\''-\^{KYCr,Y-v\^Vph{\,Y-v/2){nr{A^,Y-^/2) 

r=L,Rr]>0 

-^\(^r,Y^r,Y^r]f^ph{\,Y+r}/2)'i^r{\,Y+r}/2)^ , 
2 

{(0,F|#f |0,F + r/)) = 27ri^ ^ \Vr\''\c*^YCr,Y+^\^Vph{A^^Y+^/2) {nr{A^^Y+^/2) + l) 

r=L,R 

2 

{{0,Y\¥^^\0,Y-r])) = 27ri^ ^ |^,p|c* yC,,y_^pP^/,(A^,y_^/2)n,(A^,y_^/2) , 

^ r=L,R 

where we have used Eqs.([8| and (20). When we define the transition probabilities as 

j.Y,Y+n ^ 27r^|v^p|c*yC^,y+^pPp^(A^,y+^/2) (n^(A^,y+^/2) + 1) , 

kr'^'"^ = 27r^\Vr\'^\4^Y(^r,Y-ri\^^ph{\,Y-ri/2)nr{\,Y-r^/2) , 

where > 0, we can write 

r=L,i? 77>0 
r=L,R 

((o,r|*f |o,y-7?)) = i Y 

r=L,R 



1) 



(C9) 
(CIO) 

(Cll) 



(C12) 
(C13) 

(C14a) 
(C14b) 
(C14c) 



r 



For a TV = 2 molecule, the eigenstates of Hm are ob- eigenstates are given by 
tained as 



\1) = cl-i\L) + cji-,\R) , (C15a) 
\2) = cl-2\L) + cji-2\R) , (C15b) 



with 



= cos^ , Cji i = sin^ , (C16a) 
Cl,2 = — sin6> , Cji 2 = cos^ , (C16b) 

where tan26> = 2J/{Er — El). The eigenvalues of these 



sl^sr- +4J2 



(C17a) 



E, = ^^^^n^V(^^-^nr±±^ . (ci7b) 
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The transition probabilities are then given by 





T^9 vl 


s\v? 


26 Vph{A) (nL(A) ^ 


-1) , (C18) 


^.1,2 _ 
'^R — 




sin^ 


20 VphiA) (uRiA) - 


M) ,(C19) 


21)? 


= 


2m 


sin^ 


29 Vph{A)nL{A) , 


(C20) 


f,2,l _ 
'^R — 




sin^ 


29 Vph{A)nR{A) . 


(C21) 


2Tfi 



Appendix D: Energy Flow in terms of a energy 
change of the thermal bath 

In this section, we shall show that the formula of the 
energy flow Eg. ([61]) is also obtained by investigating the 
energy change of the thermal bath systems. We can de- 
note the energy flow from the left thermal bath as a en- 
ergy change of the left thermal bath: 

w^t) = = Y.u:L,^T^[NL,^CMBp{t)\ , (Dl) 

q 

where A/'/^q = q^L,q- This deflnition coincides with 
the ordinary deflnition of the particle current from a par- 
ticle bath to a mesoscopic system. [37] Using Eqs.(A12) 
and replacing p(t) with we have the expression of 



the energy flow from the left bath in the nonequilibrium 
stationary state as 

< = ^EE-W.--,vTr[(bt .^-bL„)Fr] . 

q ri,Y 

Substituting Eq.(49) into Eq.(D2) and after some calcu- 
lation, we obtain 

r7>0 

X |^ynL(A^^y+^/2) - ^y+r7 (^L ( A^,y+^/2 ) + l) | • 

(D3) 



The energy flow going out to the right thermal bath has 
been calculated similarly. By summing up the two con- 
tributions to obtain the same results as given in Eq.(61). 



These derivations mentioned here is to clarify that the 
energy flow which is carried by a phonon particle flowing 
from the left bath to the right one. This is the reason 
why the energy flow in the present case is described by 
the Landauer formula which is mostly used to represent 
the electronic current in which an electron flows in from 
an electron reservoir at one end and going out to the 
other end through a one-dimensional mesescopic system. 
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